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Abstract 

Deep-inelastic scattering data in the range 150 < < 35 000 GeV^ are used to investigate 
the minimum jet separation necessaiy to allow accurate description of the rate of dijet 
production using next-to-leading order perturbative QCD calculations. The required jet 
separation is found to be small, allowing about 1/3 of DIS data to be classified as dijet, 
as opposed to approximately 1/10 with more typical jet analyses. A number of precision 
measurements made using this dijet sample are well described by the calculations. The data 
are also described by the combination of leading order matrix elements and parton showers, 
as implemented in the QCD based Monte Carlo model RAPGAP. 



To be submitted to the European Physical Journal 



C. Adloff^^, V. Andreev^^, B. Andrieu^^, T. Anthonis^, V. Arkadov^^, A. Astvatsatourov'^^ 

A. Babaev23, J. Bahr'^^ P. Baranov^^, E. Barrelet^^, W. Bartel^", P. Bate2\ A. Beglarian^^, 
O. Behnke^^ C. Beieri\ A. Belousov^^, T. Benischi", Ch. Berger\ T. Bemdtl^ J.C. Bizot^^, 
V. Boudry^"^, W. Braunschweig^, V. Brisson^^, H.-B. Broker^, D.P. Brown^°, W. Briickner^^, 

D. Bruncko^^, J. Burger^°, F.W. Busser^\ A. Bunyatyan^^'^^, A. Bu^■age^^ G. Buschhorn^s, 
L. Bystritskaya^^, A.J. Campbell^°, J. Cao^^, S. Caron\ D. Clarke^ B. Clerbaux^, C. CoUard^, 
J.G. Contreras^'^\ Y.R. Coppens^ J.A. Coughlan^ M.-C. Cousinou^^, B.E. Cox^^, G. Cozzika^ 

J. Cvach29, J.B. Daintonl^ W.D. Dau^^, K. Daum'^^ •'^^ M. Davidsson^^, B. Delcourt^^, N. Delerue^^, 
R. Demirchyan3\ A. De Roeck^^ '^^ E.A. De Wolf^, C. Diaconu^^ J. Dingfelder^^ P Dixon^^ 
V. Dodonov^^ J.D. Dowell^, A. Droutskoi^^, A. Dubak^^ C. Duprel^ G. Eckerlin^°, D. Eckstein^^ 
V. Efremenko23, S. EgU^^ R. Eichler^^, F. Eisele^^ E. Eisenhandler^^, M. EUerbrock^^ E. Elsen^°, 
M. Erdmanni°'^° % W. Erdmann^^, PJ.W. Faulkner^ L. Favart^ A. Fedotov^^, R. Felst^^, J. Ferencei^°, 
S. Ferron^'', M. Fleischer^°, Y.H. Fleming'^, G. Fliigge^, A. Fomenko^"', 1. Foresti^^, J. Fomianek'^^, 
J.M. Foster2\ G. Franke^°, E. Gabathuler^^, K. Gabathuler^^, J. Garvey^ J. Gassner^^ J. Gayler^°, 
R. Gerhards^o, C. Gerlich^^ S. Ghazaryan^'^^, L. Goerlich^, N. Gogitidze^^ M. Goldberg^^, 
C. Grab^^, H. Grassier^ T. Greenshaw^^, G. Grindhammer^^, T. Hadig^^ D. Haidt^°, L. Hajduk^ 
WJ. Haynes^, B. Heinemann^^, G. Heinzelmann^^, R.C.W. Henderson^'', S. Hengstmann^^, 
H. Henschel3^ R. Heremans\ G. Herrera^''^^ I. Herynek^^, M. Hildebrandt^^ M. Hilgers'^^ 
K.H. Hiller3^ J. Hladky^g^ p noting^ D. Hoffmann^^, R. Horisberger^^ S. Hurlingi°, M. Ibbotson^^, 

i§sever^, M. Jacquet^^, M. Jaffre^^, L. Janauschek^^, X. Janssen^, V. Jemanov^^, L. Jonsson^°, 
C. Johnson^, D.P Johnson^ M.A.S. Jones^^ H. Jung^o-^o, H.K. Kastli=^^ D. Kant^^ M. Kapichine^ 
M. Karlsson^o, O. Karschnick", F Keil^"^, N. Keller^^ J. Kennedy^^ I.R. Kenyon^ S. Kermiche^^ 

C. Kiesling^s, p Kjellberg^o, M. K\em^\ C. Kleinworti°, T. Kluge\ G. Knies^^, B. Koblitz^^, 
S.D. Kolya^i, V. Korbel^^, P Kostka3^ S.K. Kotelnikov^^, R. Koutouev^^ A. Koutov®, H. Krehbiel^^, 
J. Kroseberg'^' , K. Kruger^", A. Kupper^^, T. Kuhr^^, T. Kurca^^'^^, R. Lahmann^°, D. Lamb^, 
M.PJ. Landon^^ W. Lmge^\ T. Lastovicka=^°■3^ P Laycockl^ E. Lebailly^^, A. Lebedev^^, 

B. LeiBner^, R. Lemrani^°, V. Lendermann^, S. Levonian^", M. Lindstroem^°, B. List^^, E. Lobodzinska^°'^, 
B. Lobodzinski'''^°, A. Loginov^^, N. Loktionova^"^, V. Lubimov^^, S. Liiders^^, D. Luke'^'^°, 

L. Lytkin^^ H. Mahlke-Kruger^^, N. Malden^^, E. Malinovski^"^, I. Malinovski2^ R. Maracek^^ 
P Marage^, J. Marks^^, R. MarshalF\ H.-U. Martyn\ J. Martyniak^ S.J. Maxfield^^ D. Meer3^ 
A. Mehta^^, K. Meier^^, A.B. Meyer^\ H. Meyer^^, J. Meyer^°, P-O. Meyer^, S. Mikocki^ 

D. Milstead^^, T. Mkrtchyan^^, R. Mohr^^, S. Mohrdieck^\ M.N. Mondragon^, F. Moreau^^, 
A. Morozov^, J.V. Morris^, K. MuUer'^'^, P. Murin^'''"^^, V. Nagovizin^^, B. Naroska^\ J. Naumann'', 
Th. Naumann^s, G. Nellen^^, PR. Newman^ T.C. NichoUs^ F. Niebergall^\ C. Niebuhr^°, 
O. Nix^^, G. Nowak^, J.E. 01sson^°, D. Ozerov^^ V. Panassik^, C. Pascaud^^, G.D. PateF^, 
M. Peez22, E. Perez'', J.P Phillips^^ D. Pitzl^°, R. PoschP^, I. Potachnikova^^ B. Povh^^ 
K. Rabbertz\ G. Radel\ J. Rauschenberger", P Reimer^^ B. Reisert^^ D. Reyna^^, C. Risler^^, 

E. Rizvi\ P Robmann^^ R. Roosen"^, A. Rostovtsev^^, S. Rusakov^^, K. Rybicki^ D.PC. Sankey^ 
J. Scheins\ F-P SchiUingio, P Schleperio, D. Schmidt^^, D. Schmidt^", S. Schmidt^^ S. Schmitti°, 
M. Schneider22, L. Schoeffel^ A. Schoning^e, T. Schomer^^, V. Schroder^°, H.-C. Schultz- 
Coulon^, C. Schwanenberger^°, K. Sedlak^^ F. Sefkow^^, V. Shekelyan^^, I. Sheviakov^^, 
L.N. Shtarkov^^, Y. Sirois^^, T. Sloan^^ P Smimov^^, Y. Soloviev^^, D. South^i, V. Spaskov^ 
A. Specka^^ H. Spitzeri\ R. Stamen^ B. Stella^i, J. Stiewei"^, U. Straumann^^ M. Swartl^ 

M. Tasevsky^^, V. Tchemyshov^'^, S. Tchetchelnitski^^, G. Thompson^^, PD. Thompson^, N. Tobien^", 
D. Traynor^^, P. Truol^^, G. Tsipolitis^°'^^, I. Tsurin^^, J. Tumau^, J.E. Tumey^^, E. Tzamariudaki^^, 
S. Udluft^s, M. Urban^^, A. Usik^^, S. Valkar^o, A. Valkarova^^, C. Vallee^^, p Van Mechelen^ 



1 



S. Vassmev^ Y. Vazdik^^, A. Vichnevski^ K. Wacker^ R. Wallny^^, B. Waugh^^, G. Weher^\ 
M. Weber^^, D. Wegener^ C. Wemer^^ M. Wemer^^ N. Wemer^^, G. White^^, S. Wiesand^^, 
T. Wilksen^o, M. Winders, G.-G. Winter^o, Ch. Wissing^, M. Wobisch^o, E. Wiinschio, A.C. Wyatt^i, 
J. Zacek^o, J. Zalesak^o, Z. Zhang^^, A. Zhokin^^ F. Zomer^^, J. Zsembery^ andM. zurNedden^^ 



^ /. Physikalisches Institut der RWTH, Aachen, Germany"' 
^ III. Physikalisches Institut der RWTH, Aachen, Germany^ 

^ School of Physics and Space Research, University of Birmingham, Birmingham, UK'' 

^ Inter- University Institute for High Energies ULB-VUB, Brussels; Universitaire Instelling 

Antwerpen, Wilrijk; Belgium'^ 

^ Rutherford Appleton Laboratory, Chilton, Didcot, UK'' 
^ Institute for Nuclear Physics, Cracow, Poland'^ 
^ Institut fUr Physik, Universitdt Dortmund, Dortmund, Germany'^ 
^ Joint Institute for Nuclear Research, Dubna, Russia 
9 CEA, DSM/DAPNIA, CE-Saclay, Gif-sur-Yvette, France 
DESY, Hamburg, Germany 

Institut fUr Experimentalphysik, Universitdt Hamburg, Hamburg, Germany"^ 

Max-Planck-Institut fUr Kemphysik, Heidelberg, Germany 

Physikalisches Institut, Universitdt Heidelberg, Heidelberg, Germany" 

Kirchhoff- Institut fur Physik, Universitdt Heidelberg, Heidelberg, Germany" 

Institut fUr experimentelle undAngewandte Physik, Universitdt Kiel, Kiel, Germany 

Institute of Experimental Physics, Slovak Academy of Sciences, Kosice, Slovak Republic^'^ 

School of Physics and Chemistry, University of Lancaster, Lancaster, UK'' 

Department of Physics, University of Liverpool, Liverpool, UK'' 

Queen Mary and Westfield College, London, UK'' 

Physics Department, University of Lund, Lund, Sweden^ 

Physics Department, University of Manchester, Manchester, UK'' 

22 CPPM, CNRS/IN2P3 - Universite Mediterranee, Marseille - France 

2^ Institute for Theoretical and Experimental Physics, Moscow, Russia' 

2^ Lebedev Physical Institute, Moscow, Russia'^''^ 

2^ Max-Planck-Institut fUr Physik, MUnchen, Germany 

2^ LAL, Universite de Paris-Sud, IN2P3-CNRS, Or say, France 

2^ LPNHE, Ecole Poly technique, IN2P3-CNRS, Palaiseau, France 

28 LPNHE, Universites Paris VI and VH, IN2P3-CNRS, Paris, France 

2^ Institute of Physics, Academy of Sciences of the Czech Republic, Praha, Czech Republic'^'^ 
Faculty of Mathematics and Physics, Charles University, Praha, Czech Republic'^'^ 
Dipartimento di Fisica Universitd di Roma Tre and INFN Roma 3, Roma, Italy 

^2 Paul Scherrer Institut, Villigen, Switzerland 

Fachbereich Physik, Bergische Universitdt Gesamthochschule Wuppertal, Wuppertal, Ger- 
many 

Yerevan Physics Institute, Yerevan, Armenia 
DESY, Zeuthen, Germany 

Institut fur Teilchenphysik, ETH, Zurich, Switzerland^ 
Physik-Institut der Universitdt Zurich, Zurich, Switzerland^ 



2 



Also at Physics Department, National Technical University, Zografou Campus, GR-15773 

Athens, Greece 

Also at Rechenzentrum, Bergische Universitat Gesamthochschule Wuppertal, Germany 
Also at Institutfur Experimentelle Kernphysik, Universitat Karlsruhe, Karlsruhe, Germany 
Also at Dept. Fis. Ap. CINVESTAV, Merida, Yucatan, Mexico^ 
Also at University of P.J. Safdrik, Kosice, Slovak Republic 
^'^ Also at CERN, Geneva, Switzerland 
Also at Dept. Fis. CINVESTAV, Mexico City, Mexico'' 



" Supported by the Bundesministerium fUr Bildung und Forschung, FRG, under contract num- 
bers 05 HI IGUA /I, 05 HI IPAA /I, 05 HI IPAB /9, 05 HI IPEA /6, 05 HI IVHA /7 and 05 
HI lVHB/5 

Supported by the UK Particle Physics and Astronomy Research Council, and formerly by the 
UK Science and Engineering Research Council 
" Supported by FNRS-NFWO, IISN-IIKW 

Partially Supported by the Polish State Committee for Scientific Research, grant no. 2P0310318 
and SPUB/DESY/P03/DZ-1/99, and by the German Bundesministerium fiir Bildung und Forschung, 
FRG 

^ Supported by the Deutsche Forschungsgemeinschaft 

f Supported by VEGA SR grant no. 2/1169/2001 

^ Supported by the Swedish Natural Science Research Council 

Supported by Russian Foundation for Basic Research grant no. 96-02-00019 
* Supported by the Ministry of Education of the Czech Republic under the projects INGO- 
LAll 6/2000 and LN00A006, by GA AVCR grant no Bl 01 0005 and by GAUK grant no 1 73/2000 
^ Supported by the Swiss National Science Foundation 

Supported by CONACyT 
' Partially Supported by Russian Foundation for Basic Research, grant no. 00-15-96584 



3 



1 Introduction 



One of the most remarkable results arising from the study of deep-inelastic ep scattering (DIS) 
at HERA is the large range of squared four-momentum transfer Q'^ over which perturbative 
QCD calculations are able to describe the measurements of the inclusive cross section. Next- 
to-leading order (NLO) calculations are successful from values of as low as a few GeV^ up 
to Q'^ ~ 35 000 GeV^ Investigations of the hadronic final state have shown that QCD is 

also able to describe events containing two highly energetic jets . Such investigations have 
tended to require large inter-jet separations, that is, a large relative jet transverse momentum, or 
a large transverse jet energy in the Breit frame [0,^. Typically only about a tenth of the DIS 
sample is then classified as dijet events. These large scales are chosen to avoid the region in 
which multiple parton emission is likely to become significant. Here, we examine the possibility 
that fixed order perturbative QCD is able to describe the hadronic final state in DIS even where 
jet separations are small. Some hints that this may be possible have been seen in measurements 
of event shape variables [[TO|,[TT]]. 

The study proceeds by first identifying the minimum inter-jet separation for which the rate 
of dijet production is successfully described by NLO QCD calculations. Using this separation, 
about 1/3 of DIS events are classified as dijet. Based on this sample a number of measure- 
ments is made and compared with perturbative QCD. At the lower end of the range studied, 

> 150 GeV^, the sample is dominated by gluon induced events, eg — egg, whereas at the 
upper end quark induced events predominate, eg eqg. The data thus provide a thorough test 
of the QCD calculations. 

An alternative QCD based description of the measurements is provided by Monte Carlo 
models which describe DIS using leading order (LO) QCD matrix elements matched to par- 
ton showers. The inclusion of the latter would suggest that these might describe dijet data in 
the region in which jet separations are very small. Indeed, studies have been made of sub-jet 
multiplicities and jet shapes, in which Monte Carlo models incorporating parton showers have 



performed reasonably well [|12|,|130. The LO calculations incorporated in these models should 
ensure they are also able to describe data at large jet separations. However, their overall per- 
formance in describing the hadronic final state in DIS has not yet been satisfactory Q]. In 
the present analysis, we confront our measurements with the QCD model RAPGAP [[T^] which 
was not considered in [|T^. 



2 Experimental procedure 
2.1 Selection of DIS events 

The present analysis is based on a data sample corresponding to an integrated luminosity of 
~ 35 pb^^ recorded in 1995-97 with the HI detector at HERA. In this period HERA operated 
with positron and proton beams of 27.5 GeV and 820 GeV energy, respectively, yielding a 
centre-of-mass energy ^/s of 300 GeV. 

A detailed description of the HI detector is given in [|1^. The detector components of 
most importance for this study are the central tracking system and the liquid argon calorimeter. 
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We use a coordinate system with its origin at the nominal interaction point and its positive z 
axis along the direction of the outgoing proton beam. Polar angles are denoted by Q and the 
"forward" region is that with Q < 90°. Neutral current DIS events are selected using criteria 
similar to those described in [^. These include the requirement that a scattered positron be 
identified in the liquid argon calorimeter at a polar angle 9e < 150°. The positron reconstruction 
method and the fiducial cuts on the positron impact position in the liquid argon calorimeter are 
described in The value of Q^, determined from the energy and polar angle of the scattered 
positron, must exceed 150 GeV^. The Bjorken scaling variable y must satisfy 0.1 < y < 0.7. 
At low y, measurements of y using the polar angles of the positron and of the reconstructed 
hadronic final state [|l]], y^a, are more precise than those using the energy and polar angle of the 
scattered positron, y^,. At high y the situation is reversed. The low y requirement is thus applied 
using the double angle measurement, 0.1 < yda, whereas the high y restriction is applied using 
the positron measurement, y^. < 0.7. The selection yields a sample of ~ 60 000 DIS events with 
negligible background [0] . 



2.2 Jet algorithm and observables 

Jets are reconstructed with the modified Durham algorithm, described in more detail in [^|T7]], 
which is applied in the laboratory frame. Hadronic energy deposits measured in the liquid argon 
calorimeter and the backward "SpaCal" calorimeter are used, as are tracks reconstructed in the 
central tracking chambers, avoiding double counting of energy. All of these are referred to as 
"proto-jets" in the following and are required to have a polar angle greater than 7° to ensure 
that they are well measured. The proton remnant, which escapes direct detection, is included 
in the jet reconstruction by forming a missing-momentum four-vector, which is treated as an 
additional proto-jet. 

The algorithm uses the relative k'^^j = 2rmn[Ef,Ej] (1 — cos%) of proto-jets i, j as a 
measure of their separation, where Ei and Ej are the energies of the proto-jets i and j, and 
9ij the angle between them. The pair i, j with the minimum kxij is combined to form a new 
proto-jet by adding the four-momenta pi and pj. The iterative clustering procedure is ended 
when exactly two final state jets and the proton remnant jet remain. 

In order to select a sample of dijet events we define the variable y2: 



where i, j may be any of the two final state jets or the remnant jet, and W is the invariant mass 
of all objects entering the jet algorithm, including the missing-momentum vector. It is ensured 
that the jets are well contained within the liquid argon calorimeter by requiring that, for both 
non-remnant jets, 10° < 6'jet < 140°. 

We study the following observables: ?/2 as defined above; the polar angles of the forward 
and the backward (non-remnant) jets in the laboratory frame 6'fwd and 9bv,d\ the dimensionless 
variables Xp and Zp-, and the average transverse energy of the two final state (non-remnant) jets 
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in the Breit frame i?TBrcit- In order to calculate ExBmiu the jets found in the laboratory frame 
are boosted into the Breit frameQ. We calculate Zp and Xp according to 

^ mmi=i^2Ei (1 - cos 9i) ^ 

where Ei and 9i are the energies and polar angles of the two (non-remnant) jets, and mi2 is the 
invariant dijet mass. Matrix elements for dijet production are frequently expressed using these 
variables [[r^-|20|]. In leading order QCD, the cross section diverges for — > and Xp — > 1 due 
to coUinear and infrared singularities. 

In the QCD models and the NLO calculations, jets are defined by applying the above algo- 
rithm to the four- momenta of hadrons or partons. In particular, the requirement that the polar 
angle be greater than 7° is always applied. 



2.3 Data correction and systematic uncertainties 

The procedures used for data correction and the determination of the systematic uncertainties 
are similar to those described in 0]. The measured jet distributions are corrected for the 
effects of the limited detector acceptance and resolution, and for the effects of QED radiation. 
This is done using bin-by-bin correction factors determined with the QCD Monte Carlo mod- 
els ARIADNE [P] and LEPTO [|2|], both of which are incorporated in DJANGO The 



average of the measured jet distributions corrected with ARIADNE or LEPTO is taken as the 
final result. As a cross check, some distributions are also corrected using a regularized unfold- 



ing technique [ p4| ] (for a brief description see [pSQ). The two correction methods lead to very 
similar results. 

The dominant systematic errors are due to the model dependence of the corrections and the 
uncertainty of the electromagnetic and hadronic energy scales of the calorimeters. The errors 
are added in quadrature to yield the total systematic error. The model uncertainty is taken to 
be the difference between the averaged correction factors and those determined with a single 
model, and is of the order of 5%. The uncertainty of the electromagnetic energy scale of the 
liquid argon calorimeter ranges from ±0.7% to ±3%, depending on the scattered positron's 
impact position. The changes in the measured dijet distributions resulting from the variation of 
this scale within its uncertainty lead to a systematic error of less than 1%. The hadronic energy 
scale of the liquid argon calorimeter is varied by ±4%, which leads to an average uncertainty 
of 4% in the dijet measurements. 



2.4 Perturbative QCD and model calculations 

The perturbative QCD predictions presented in the following are calculated using the DISENT 



program [pq]. The agreement of DISENT with other NLO programs is discussed in [p7|-p0|]. In 



'The Breit frame is related to the hadronic centre-of-mass frame by a longitudinal boost. In both frames, the 
total transverse momentum of the hadronic final state is zero whereas in the laboratory frame it is constrained to 
balance the scattered lepton's transverse momentum. 
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the calculations, we use the CTEQ5M parton density functions [Pl]], choose Q as the renormal- 
ization and factorization scale if not otherwise stated and set the value of as{Mz) to 0.1183. 
This gives a good description of the inclusive DIS cross section in the kinematic range of this 
analysis 0]. Other choices of recent parton density parameterizations, for example those de- 
termined by the HI collaboration in are found to yield very similar NLO predictions. The 
size of the hadronization effects is determined using the QCD Monte Carlo models LEPTO and 
ARIADNE. Hadronization correction factors are obtained by dividing the jet distributions for 
hadrons by those determined from the partons after the parton shower or dipole cascade, respec- 
tively. The average of the hadronization corrections obtained from the two models is applied to 
the NLO calculations in the comparisons. The uncertainty of this procedure is conservatively 
estimated to be half the size of the correction. This is significantly larger than the difference 
between the corrections determined with ARIADNE and LEPTO. 

Comparisons are also made with the QCD based Monte Carlo program RAPGAP. This mod- 



els QCD radiation with initial and final state parton showers [ p2[ ] combined with leading order 



QCD matrix elements []TE[-PD|]. Hadronization is simulated using the Lund string model [R^IB^]. 



We use the default model parametersFI and the CTEQ4L parton density functions [BSl]. 



3 Determination of minimum required jet separation 

A direct way of determining the minimum jet separation necessary to ensure that NLO QCD 
describes the dijet production rate is to compare measurements of the jet separation itself with 
calculations. The measured ?/2 distribution, normalized to the inclusive DIS cross section 
in the region defined by Q"^ > 150 GeV^, 0.1 < y < 0.7 and 9^ < 150°, is shown with the 
results of various calculations in Figure [T| and listed in Table [T|. 

We observe that the NLO perturbative QCD calculations combined with hadronization cor- 
rections overestimate the measured cross section drastically in the region y2 < 0.001, where jet 
separations are smallest. Here the difference between LO and NLO predictions is large. The 
renormalization scale dependence, estimated by varying fiR in the range Q/2 < < 2Q, 
and the hadronization corrections are also large. All three criteria suggest that fixed order per- 
turbative QCD predictions are not reliable in the region y2 < 0.001, and agreement of the 
calculations with the data cannot be expected. The situation is much improved at y2 > 0.001, 
and a good description of the data is observed from y2 ~ 0.001 up to the largest y2 values, 
where jet structures are most distinct. (The deviation of NLO QCD in the highest 7/2 bin can 
be explained by the exceptionally large parton density function dependence in this region of the 
dijet phase space.) 

The 1/2 distribution of Figure |T] is also compared to the QCD model RAPGAP, which de- 
scribes the y2 cross section over the full measured range. In particular, the region of very low 
y2 is well described, which suggests that the combination of parton showers and Lund string 
hadronization used in RAPGAP accurately models multi-parton emissions. 

^The cut-off parameter for the LO matrix element calculation PT2CU is set to 5 GeV^. 
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4 Study of dijet sample 



Motivated by the agreement of NLO perturbative QCD with the data at y2 > 0.001 in Figure [T|, 
we investigate this sample of dijet events in more detail. This sample contains about 1/3 of the 
selected DIS events. 

The dijet cross section da2 / dEx Breit in several ranges of is shown in Figure |^ and 
listed in Table ^ The normalization is again to a^is, the DIS cross section in the region 
defined by 0.1 < y < 0.7, < 150° and the indicated Q"^ range. For a sizable fraction of the 
events, i^^Breit is smaller than 5GeV and the mean value of -ETBreit over the dijet sample is 
~ 6 GeV. The more restrictive dijet samples used in other QCD analyses typically require that 
one of the jets have Et Breit > 7.5 GeV or higher [^H]. We compare the measurements with 
perturbative QCD calculations in NLO for two choices of renormalization scale, fiR = Q and 
fiji = i^TBreit-PI Pcrturbative QCD in NLO describes the Etbtch distributions well, including 
the region E^Breit < 5 GeV. Although the numerical values of Q'^ and i^j^B^cit different 
for most events, the difference between the NLO predictions for the different scales is small. 
RAPGAP also describes the Et Bmit distributions well. 

The measured distributions of the forward and backward jet polar angles, 6'fwd and 6'bwd, are 
shown in Figure |^ (and Table ^ in five ranges. The ^^f^d distributions increase strongly 
towards small angles and are less dependent on Q'^ than the 6'bwd distributions, as expected 
if the forward jets are largely due to initial state radiation off the constituents of the proton. 
The distributions are well described by NLO perturbative QCD and by the RAPGAP model. 
In particular at small jet polar angles and at relatively small Q^, the LO calculations differ 
considerably from the NLO ones and are unable to describe the data. A prediction including 
only the phase space contribution, without the QCD matrix elements, results in relatively flat 
^^fwd distributions and does not describe the data. 

The 6'bwd distribution has its maximum at large polar angles in the kinematic region 150 < 
< 275 GeV^. With increasing Q"^, the maximum shifts into the forward region of the 
detector. Again, this is as expected due to the increasing fraction of the proton's momentum 
transferred to the hadronic final state as grows. Perturbative QCD in NLO and the QCD 
model RAPGAP describe the distributions well, while the LO QCD predictions fail. 

The Xp and Zp distributions are shown in Figure ^ and listed in Table ^. The Zp distribution is 
well described by both NLO QCD and RAPGAP. The NLO calculations without hadronization 
corrections are also shown. They describe the data fairly well since the hadronization correc- 
tions are small. (Note that, by definition, the error band includes half of the hadronization 
corrections.) The Xp distribution is also well described by NLO QCD. The NLO predictions 
including only the quark contribution (eg eqg) to the dijet cross section are shown sepa- 
rately. The proportion of the dijet cross section that is quark-induced is expected to be largest at 
large Xp values, and varies from ^ 30% at the lowest range to nearly 100% at > 10 000 
GeV^. This illustrates that our measurements are sensitive tests of both the gluon- and the 
quark-initiated NLO matrix elements. 

''For the choice fiji = EtbicH, the cut £'TBioit > 0.5 GeV is applied to improve the convergence of the 
DISENT calculations. This has a negligible effect in the selected dijet phase space. 
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5 Summary 



We have investigated the minimum inter-jet separation necessary to ensure that next-to-leading 
order QCD calculations are able to accurately describe dijet production in deep-inelastic scat- 
tering. Using data in the kinematic range 150 < < 35 000 GeV^, the required separation is 
found to be small, resulting in the selection of a dijet sample containing about 1/3 of the DIS 
events, a significantly larger proportion than the approximately 1/10 obtained with more typical 
jet selection criteria. Measurements of the distribution of variables sensitive to the dynamics 
of jet production are well described by NLO QCD calculations, for either choice of renormal- 
ization scale Q or ExBrdt- This good description extends to unexpectedly small jet separations 
and covers regions in which both gluon and quark induced processes dominate. 

Due to their precision and to the large phase space covered, the measurements also signif- 
icantly constrain QCD Monte Carlo models. A good description of the data may be achieved 
by models which combine leading order QCD matrix elements with parton showers, as demon- 
strated here in the case of RAPGAP. 
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1 /~\rr 11/^1 


1 da 


Sstat(%) 


Ssys(%) 




150 < g2 < 5000 GeV^ 


5.5 - 4.5 


0.00912 


±5.2 


±8.8 


4.5-3.75 


0.255 


±1.5 


±4.1 


3.75 - 3.0 


0.465 


±1.1 


±0.9 


3.0-2.45 


0.335 


±1.3 


±4.6 


2.45 - 2.3 


0.241 


±2.5 


±4.3 


2.3 - 2.0 


0.135 


±2.4 


±5.4 


2.0 - 1.7 


0.0411 


±4.1 


±8.7 


1.7-1.1 


0.0049 


±8 


±11 



Table 1: Normalized jet cross sections as a function of y2 for > 150 GeV^, 9e < 150°, 
0.1 < y < 0.7, and 10° < ^jet < 140°. The relative statistical errors 5stat and systematic errors 
6sys are given in per cent. 



Et Brcit 

[GeV] 


1 da2 


Sstat(%) 




1 da2 


Sstat(%) 




^DIS Brcit 

[GeV-1] 


""DIS d^r Brcit 

[GeV-1] 




150 < g2 < 275 GeV2 


275 < g2 < 575 GeV2 


0-5 


0.013 


±2.5 


±9.9 


0.0165 


±2.9 


±8.8 


5-10 


0.0233 


±2.2 


±5.9 


0.027 


±3 


±21 


10-17.5 


0.00619 


±3.3 


±7.4 


0.00851 


±3.6 


±9.7 


17.5-25 


0.00089 


±8 


±21 


0.0015 


±8 


±15 


25-35 


0.00012 


±18 


±17 


0.00023 


±17 


±23 




575 < g2 < 5000 GeV^ 


5000 < g2 < 10 000 GeV^ 


0-5 


0.0193 


±4 


±10 


0.035 


±29 


±39 


5-10 


0.0301 


±3.3 


±4.6 


0.029 


±23 


±27 


10 - 17.5 


0.0110 


±4.4 


±7.3 


0.0089 


±30 


±15 


17.5-25 


0.0031 


±8 


±15 


0.0040 


±42 


±14 


25-35 


0.00053 


±14 


±14 


0.0017 


±57 


±26 




Q'^ > 10 000 GeV^ 




0-5 


0.049 


±76 


±100 


5-10 


0.032 


±63 


±44 


10 - 17.5 


0.0061 


±89 


±100 


17.5-25 


0.0032 


±100 


±23 


25-35 


0.0025 


±100 


±19 



Table 2: Normalized dijet event cross sections as a function of .ETBreit, determined with the 
modified Durham algorithm. The selection criteria are 9e < 150°, 0.1 < y < 0.7, y2 > 0.001 
and 10° < 6'jet < 140°. The range is given in the table. The relative statistical errors 6stat 
and systematic errors Sgys are given in per cent. 
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ffwd 

[deg] 


1 da2 


Sstat(%) 


5sys(%) 


1 da2 


Sstat(%) 


Ssysi%) 


cdis d6'fwd 
[deg-i] 


o'Dis dOf^d 
[deg-i] 




150 < g2 < 275 GeV^ 


275 < < 575 GeV^ 


10-20 


0.00555 


±2.8 


±2.6 


0.00858 


±2.9 


±2.5 


20-35 


0.00422 


±2.6 


±2.5 


0.00592 


±2.9 


±3.8 


35 - 60 


0.00236 


±2.9 


±7.9 


0.00299 


±3.3 


±3.3 


60 - 90 


0.00129 


±3.7 


±5.8 


0.00105 


±4.9 


±4.5 


90 - 140 


0.000364 


±5.1 


±4.6 


0.00231 


±8.1 


±6.2 




575 < < 5000 GeV- 


5000 < Q- < lOOOOGeV- 


10-20 


0.0142 


±3.4 


±5.5 


0.020 


±22 


±4 


20 - 35 


0.00745 


±3.7 


±3.9 


0.0085 


±25 


±3 


35-60 


0.00292 


±4.7 


±2.7 


0.0034 


±31 


±6 


60-90 


0.000921 


±7.4 


±4.8 


0.00017 


±100 


±4 


90 - 140 


0.000146 


±15 


±8.5 





< 0.0001 (68% CL) 




g2 > lOOOOGeV^ 




10 - 20 


0.026 


±54 


±8 


20 - 35 


0.0087 


±69 


±3 


35-60 





< 0.002 


:68% CL) 


/J 

fbwd 

[deg] 


1 Clo-2 

o-Dis dObwd 
[deg-1] 


Sstat(%) 




1 do-2 
[deg-i] 


Sstat(%) 


Ssys(%) 




150 < < 275 GeV^ 


275 < < 575 GeV^ 


10-40 


0.000711 


±3.7 


±3.9 


0.00171 


±3.2 


±5.1 


40-60 


0.00203 


±3.2 


±4.6 


0.00343 


±3.3 


±3.2 


60-80 


0.00202 


±3.3 


±7.2 


0.00254 


±3.8 


±3.7 


80 - 100 


0.00221 


±3.3 


±3.7 


0.00230 


±4.0 


±3.0 


100 - 120 


0.00219 


±3.3 


±5.2 


0.00200 


±4.4 


±4.3 


120 - 140 


0.00218 


±3.3 


±4.6 


0.00171 


±4.6 


±3.5 




575 <Q^ < 5000 GeV^ 


5000 < < 10 000 GeV^ 


10-40 


0.00462 


±3.3 


±8.6 


0.0085 


±20 


±3 


40-60 


0.00387 


±4.4 


±2.6 


0.0042 


±30 


±6 


60 - 80 


0.00245 


±5.4 


±2.6 


0.0015 


±48 


±16 


80 - 100 


0.00206 


±6.0 


±3.6 


0.0011 


±52 


±11 


100 - 120 


0.00162 


±6.8 


±3.0 


0.00035 


±78 


±4 


120 - 140 


0.000959 


±8.0 


±3.8 


0.00038 


±57 


±19 




> lOOOOGeV^ 




10 - 10 


0.0099 


±52 


±6 


40-60 


0.0017 


±100 


±4 


60-80 


0.0038 


±89 


±6 


80 - 100 





< 0.002 


;68% CL) 



Table 3: Normalized dijet event cross sections as a function of 6'fwd and 6'bwd- The selection 
criteria are given in table 2. The relative statistical errors Sgtat and systematic errors 6sys are 
given in per cent. 
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1 do-2 

(^DIS dZp 


Sstat(%) 




1 do-2 

(^DIS dZp 


Sstat(%) 


Ssys(%) 




150 < < 275 GeV^ 


275 < g2 < 575 GeV^ 


- 0.125 


0.377 


±3.0 


±5.4 


0.527 


±3.4 


±6.1 


0.125 — 0.25 


0.620 


±2.5 


±5.8 


0.728 


±2.9 


±3.8 


0.25 — 0.375 


0.474 


±2.8 


±4.4 


0.576 


±3.2 


±2.2 


0.375 — 0.5 


0.412 


±3.0 


±5.2 


0.514 


±3.4 


±4.0 




575 < g2 < 5000 GeV^ 


5000 < < 10 000 GeV^ 


— 0.125 


0.681 


±4.2 


±4.1 


0.65 


±30 


±7 


0.125 - 0.25 


0.804 


±4.0 


±3.0 


0.74 


±29 


±10 


0.25 - 0.375 


0.710 


±4.1 


±7.1 


1.0 


±26 


±3 


0.375 - 0.5 


0.685 


±4.2 


±4.9 


0.88 


±27 


±7 




> 10 000 GeV^ 




0-0.125 


1.2 


±63 


±3 


0.125 - 0.25 


0.52 


±89 


±6 


0.25 - 0.375 


0.27 


±100 


±5 


0.375 - 0.5 


1.1 


±69 


±16 




1 do-2 
cdis dxp 


Sstat(%) 


Ssys(%) 


1 do-2 

(^DIS dXp 


Sstat(%) 


Ssys(%) 




150 <Q'' < 275 GeV^ 


275 <Q'^ < 575 GeV^ 


- 0.2 


0.22 


±3 


±10 


0.12 


±6 


±13 


0.2 — 0.4 


0.368 


±2.7 


±8.6 


0.338 


±3.5 


±8.2 


0.4 - 0.6 


0.317 


±2.7 


±7.3 


0.424 


±3.2 


±5.5 


0.6 — 0.8 


0.23 


±3 


±12 


0.415 


±3.0 


±6.3 


0.8 — 1.0 


0.04 


±5 


±12 


0.17 


±4 


±10 




575 < g2 < 5000 GeV^ 


5000 < g2 < lOOOOGeV^ 


- 0.2 


0.039 


±13 


±22 




0.2 — 0.4 


0.190 


±6.2 


±9.5 





< 0.03 (68% CL) 


0.4 - 0.6 


0.348 


±4.8 


±5.1 


0.11 


±52 


±29 


0.6 - 0.8 


0.573 


±3.8 


±6.4 


0.25 


±37 


±14 


0.8-1.0 


0.65 


±3 


±10 


1.7 


±18 


±4 




Q'' > lOOOOGeV^ 




0.6 - 0.8 





< 0.2 (68'X CL) 


0.8-1.0 


2.0 


±47 


±46 



Table 4: Nomialized dijet event cross sections as a function of Zp and Xp, detennined with the 
modified Durham algorithm. The selection criteria are given in table 2. The relative statistical 
errors Sgtat and systematic errors Sgys are given in per cent. 
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Figure 1: (a) Distribution of ?/2 for > 150 GeV^, determined with the modified Durham al- 
gorithm. The events satisfy 10° < ^^jet < 140°. Here, and in the following figures, the statistical 
errors are given by the inner error bars and the outer error bars correspond to the quadratic sum 
of the statistical and systematic errors. Also shown are perturbative QCD calculations in LO, in 
NLO with and without hadronization corrections, and the predictions of the QCD model RAP- 
GAP. The shaded band shows the renormalization scale uncertainty of the NLO calculations, 
which is estimated here and below by varying fin in the range Q/2to2Q. 
(b) The ratios of the data and various predictions. The vertical error bars correspond to the 
uncertainty of the data only. 
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Figure 2: (a) Dijet event Et Breit distributions determined with the modified Durham algorithm 
in various ranges. The dijet events satisfy y2 > 0.001 and 10° < ^^jct < 140°. Also 
shown are perturbative QCD calculations in NLO with = Q and fiR = ETBrdt, and the 
predictions of the QCD model RAPGAP. The shaded band corresponds to the quadratic sum of 
the hadronization and renormalization scale uncertainties. 

(b) The ratios of the data and various predictions. The vertical error bars correspond to the 
uncertainty of the data only. 
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Figure 3: Dijet event 6'twd and ^^bwd distributions determined with the modified Durham algorithm 
in various Q'^ ranges. The dijet events satisfy y2 > 0.001 and 10° < ^^jot < 140°. Also shown are 
perturbative QCD calculations in NLO and LO, the predictions of the QCD model RAPGAP, 
and a phase space calculation in arbitrary normalization. The shaded band corresponds to the 
quadratic sum of the hadronization and renormalization scale uncertainties. 
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Figure 4: Dijet event Zp and Xp distributions determined with the modified Durham algorithm 
in various ranges. The dijet events satisfy y2 > 0.001 and 10° < 6'jct < 140°. Also 
shown are perturbative QCD calculations in NLO with and without hadronization corrections, 
and the predictions of the QCD model RAPGAP. The shaded band corresponds to the quadratic 
sum of the hadronization and renormalization scale uncertainties. Note that the quark-induced 
contribution to the dijet cross section is close to 100% for Q'^ > 5000 GeV^. 
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